Testing context uality on quantum ensembles with one clean qubit. 



O. Moussa,i'0C. A. Ryan/ D. G. Cory,^'^ and R. Laflamme^' ^ 

^Institute for Quantum Computing and Department of Physics and Astronomy, 
University of Waterloo, Waterloo, ON, N2L3G1, Canada. 
^ Department of Nuclear Science and Engineering, MIT, Cambridge, Massachusetts 02139, USA 
^ Perimeter Institute for Theoretical Physics, Waterloo, ON, N2J2W9, Canada 

We present a protocol to evaluate the expectation value of the correlations of measurement out- 
comes for ensembles of quantum systems, and use it to experimentally demonstrate-under an as- 
sumption of fair sampling-ihe violation of an inequality that is satisfied by any non-contextual 
hidden-variables (NCHV) theory The experiment is performed on an ensemble of molecular nu- 
clear spins in the solid state, using established Nuclear Magnetic Resonance (NMR) techniques for 
quantum information processing (QIP). 



The Bell-Kochen-Specker (BKS) theorem [l El HIl] 
states that no noncontextual hidden-variables (NCHV) 
theory can reproduce the predictions of quantum me- 
chanics for correlations between measurement outcomes 
of some sets of observables. Any such set of observables 
constitutes a proof of the theorem. Recently, Cabello [5^ 
and others [6 used BKS proofs to derive a set of in- 
equalities that are satisfied by any NCHV theory but are 
violated by quantum mechanics for any quantum state. 
These inequalities bound certain linear combinations of 
ensemble averages of correlations between measurement 
outcomes of compatible observables; thus creating a sep- 
aration between the predicted outcomes of quantum me- 
chanics, and the bound that is satisfied by NCHV theo- 
ries. 

This provides an opportunity to test noncontextual- 
ity with finite-precision experiments -which has been 
the subject of contention for many years [71 El [9]- and 
without the need for the creation of special quantum 
states [ini HH US] • Already, two experiments, on a pair 
of trapped ^^Ca+ ions [TS , and with single photons [M], 
have demonstrated this state-independent conflict with 
noncontextuality. In this letter, we examine testing con- 
textuality on quantum ensembles. 

This manuscript is organized as follows. First, we 
sketch the arguments leading to one of the inequalities 
derived in [5^ . Then we present an algorithm to estimate 
the expectation value of the correlations of measurement 
outcomes for ensembles of quantum systems. And lastly, 
we report and discuss the result of experimentally imple- 
menting the algorithm on a 3-qubit ensemble of molecular 
nuclear spins in the solid state. 

Inequality — For a quantum system prepared ac- 
cording to some state, p, one can assign simultaneous 
outcomes {^{Sk)} of measurements of a set {Sk} of 
coobservables (i.e. comeasureable; mutually compatible; 
commuting). In this case, the correlation between the 
measurement outcomes is given by 

k k 

irrespective of the product ordering. Repeating the 



preparation and measurement many times, and averaging 
over the outcomes, one obtains an estimate of the ensem- 
ble average of the correlation {7T^Sk})p — (life ^{^k)) p • 

For the case where the coobservables {Sk} are also di- 
chotomic, with possible outcomes {v^Sk) = ±1}, the cor- 
relation ([T]) also takes on the possible values ±1, and the 
ensemble average satisfies — 1 < {^{Sk})p ^ +1- Note, 
that in this case, these operators are Hermitian and uni- 
tary (also known as Quantum Boolean Functions). 

Consider any set of observables with possible outcomes 
±1 arranged in a 3 x 3 table such that the observables in 
each column and each row are coobservable. It has been 
shown [5 that, for any NCHV theory, 

P=M^{^r,)^M^{^c,)^{^c,)-{^c,) <4, (2) 

where (tt^i) is the ensemble average of the correlation 
between outcomes of the observables listed in the first 
row, and so forth. The above inequality is independent 
of the preparation of the ensemble, provided all terms are 
estimated for the same preparation. 
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TABLE I: List of the 2-qubit observables used to show that 
quantum mechanics violates inequality ([2| . This list has been 
used by Peres [15] and Mermin [16] as a BKS proof for 4 
dimensional systems. {1,X,Z,Y} are the single qubit pauli 
operators, and, e.g. ZX := Z X indicates a measurement 
of the pauli-Z on the first qubit and pauli-X on the second. 

Now, consider a 2-qubit system (e.g. 2 spin-^ parti- 
cles), and the set of observables listed in Table |lj For 
any NCHV theory, the inequality ([2| holds for the cor- 
relations between measurement outcomes of the coob- 
servables listed in each row and column, where, for e.g., 
(tt^,) = (TTj^i^i^,^^}) = {7L1 ■ 17L • 7L7L), and so forth. 
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On the other hand, according to quantum mechanics, 
the ensemble average {7T{Sk})p given by tr(p 11^^^)- 
Thus, for a set of coobservables whose product is pro- 
portional to the unit operator -as is the case for all rows 
and columns of Table |l]- the quantum mechanical predic- 
tion of the ensemble average of the correlation is equal 
to the proportionality constant, independent of the initial 
preparation of the system. Hence, the quantum mechan- 
ical prediction for p is 6, which violates inequality 

Algorithm — To measure the correlation between 
a set of coobservables, consider introducing an ancillary 
(probe) qubit, and applying a transformation to the 
composite system for each observable Sk^ in a manner 
reminiscent of coherent syndrome measurement in quan- 
tum error correction [17 . For an observable S with the 
spectral decomposition S = P-^ — P-, where P+ and P_ 
are the projectors on the +1 and —1 eigenspaces of 5, the 
transformation L^^ is defined as Us = I2 ^ P-\- + ^ (8) P- . 
That is to say, if the system is in a -1 eigenstate of S, ap- 
ply a phase fiip (pauli-Z) to the probe qubit, and if it is 
in a +1 eigenstate, do nothing. This transformation can 
also be expressed as a controlled operation dependent on 
the state of the probe qubit: 

[/5 = I2 P+ + ^ P- 

= ^ (I2 + (P+ + P-) + ^ (I2 - (P+ - P-) 
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which is unitary for S unitary. If the system is initially 
prepared according to p, and the probe qubit in the +1 
eigenstate of the pauli-X operator, |+) = (|0) + |l))/>/2, 
the possible outcomes of pauli-X measurement on the 
probe qubit is ±1, with probabilities given by: 

p{±l) =tra+s [Us{\+){+\<i^p)Ul (|±X±| ® Id)" 

=tr, [(±|(l2 ® P+ + Z ® P-)(|+X+| ® I ±)] 

=tr« [P±p ] , 

and the ensemble average 



=+p(+l)-p(-l) 

= tr, [P+p]-tr, [P.p] 



(3) 



Thus, to measure the ensemble average of the correla- 
tion between a set of coobservables, one prepares a probe 
qubit in the +1 eigenstate of X, and the system accord- 
ing to p. As shown in Figure [l] one then applies the 
unitaries Sk in succession to the system, controlled on 
the state of the probe qubit. Since, by definition, all Sk 
mutually commute, then the order of their application 
has no bearing on the measurement outcome. Repeating 
this procedure, and averaging the outcome of the mea- 
surement on the probe system produces the correlation 
between this set of observables. Alternatively, one could 
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FIG. 1: A quantum network to encode the correlation be- 
tween the outcomes of measurements {Sk}k^i...m on a d- 
dimensional system, in the phase of a probe qubit state. Re- 
peating this procedure for the same preparation p and aver- 
aging the outcome of the measurement on the probe qubit 
gives the ensemble average {S1S2 ■ ■ ■ Sm) p- Alternatively, for 
an ensemble of quantum systems initially prepared according 
to p, on which operations are applied in parallel to the in- 
dividual systems, an ensemble measurement readily produces 

{S1S2 • • • Sm) p- 



prepare an ensemble of systems according to p; apply 
the transformations [/^ in parallel to each member of the 
ensemble; and perform a bulk ensemble measurement to 
estimate {'T^{Sk})p- This alleviates the need for isolation 
of single quantum systems, and the repeated application 
of single shot, projective measurement. 

Since inequality ([2| is valid for any preparation p, then 
one is free to choose to prepare the system according 
to the maximally mixed state. In which case, only one 
qubit -the probe system- is not maximally mixed. This 
corresponds to the model of computation known as De- 
terministic Quantum Computation with one clean qubit 
(DQCl) [18 . 

Two models — Suppose the measurement process on 
the probe qubit was e-efficient, i.e. returning a faith- 
ful answer e fraction of the time, and otherwise a uni- 
formly distributed random outcome. The probabilities 
p(±l) of obtaining outcomes ±1 will be modified to: 
= + e tr^ [P±p ], and the ensemble average 
to (X Id) = e {S)p. One can then estimate the expec- 
tation value {S)p under an assumption of fair sampling 
and knowing the value of e, which can be established from 
e(l)p. This model is equivalent to one where the probe 
system is initially in the mixed state (1 — e)^ e|+)(+|, 
provided the reduced dynamics on the probe qubit from 
preparation to measurement is represented by a unital 
map; i.e. a map that preserves the totally mixed state. 
To see this, suppose we prepare the probe qubit in some 
state pa and then apply some transformation to the com- 
posite system, whose reduced dynamics on the probe 
qubit is described by a unital linear map A. An e-efficient 
measurement of X = |+)(+| — |— )(— | has two possible 
outcomes ±1 with probabilities 
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FIG. 2: The quantum networks for the six experiments to 
estimate /3 as given in The ensemble is initially prepared 
according to pe ^ (g) where pe — (1 — e)^ + e|+X+|, 
and ^ is the single-qubit maximally mixed state. 
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FIG. 3: Malonic acid (C3H4O4) molecule and Hamiltonian 
parameters (all values in kHz). Elements along the diagonal 
represent chemical shifts, 0;^, with respect to the transmitter 
frequency (with the Hamiltonian ^-TiijOiWji). Above the di- 
agonal are dipolar coupling constants (5^^<^^ vrDi,^ (2 — 
ILilLj — YiYj), and below the diagonal are J coupling con- 
stants, {J2^<J f J*,j(^i^j +XiX^- + YiY^-). An accurate nat- 
ural Hamiltonian is necessary for high fidelity control and is 
obtained from precise spectral fitting of (also shown) a proton- 
decoupled ^^C spectrum following polarization-transfer from 
the abundant protons. The central peak in each quintuplet is 
due to natural abundance ^^C nuclei present in the crystal at 
~ 1%. (for more details see J^i20j and references therein.) 



which are precisely the statistics one obtains in case the 
probe qubit is initially in the state {1 — e)^ -\- e pa^ and 
the measurement process is faithful. 

Experiment — We implement the algorithm described 
above to perform an experimental measurement of the 
correlations as described in inequality Q on an ensemble 
of nuclear spins in the solid state using established NMR 
techniques for QIP \W,'20' . Figure |2] shows the SIX ex- 
periments required to estimate the six terms in The 
pulse sequence implementing the measurement of some 
observable is the same whether it is being measured with 
the coobservables listed in its row or column. 

The experiments were performed in a static field of 
7. IT using a purpose-built probe. The sample is a macro- 
scopic single crystal of Malonic Acid (C3H4O4), where a 
small fraction (~ 3%) of the molecules are triply labeled 
with ^^C to form an ensemble of processor molecules. 
During computation, these processors are decoupled from 
the 100% abundant protons in the crystal by applying a 
decoupling pulse sequence [21j to the protons. Shown 
in Figure [3] is a proton-decoupled ^^C spectrum, follow- 
ing polarization-transfer from the abundant protons, for 
the particular orientation of the crystal used in this ex- 
periment. A precise spectral fit gives the Hamiltonian 
parameters (listed in the inset table in Figure |3|, as well 
as the free-induction dephasing times, T2*, for the vari- 
ous transitions; these average at ~ 2ms. The dominant 



contribution [19] to T2 is Zeeman-shift dispersion, which 
is largely refocused by the control pulses. Other contri- 
butions are from intermolecular ^^C-^^C dipolar coupling 
and, particularly for Cm, residual interaction with neigh- 
boring protons. 

The carbon control pulses are numerically optimized 
to implement the required unitary gates using the 
GRAPE |22j algorithm. Each pulse is 1.5ms long, and is 
designed ^ to have an average Hilbert-Shmidt fidelity 
of 99.8% over appropriate distributions of Zeeman-shift 
dispersion and control-fields inhomogeneity. 

The two spin-^ nuclei Ci and C2, constituting the sys- 
tem on which the measurements are performed, are ini- 
tially prepared according to the totally mixed state. C^, 
representing the probe qubit, is initially prepared accord- 
ingtop, = + where |+) = ^(|0) + |1)), 

and the computational basis, {|0), |1)}, is the eigenbasis 
of the single spin Zeeman Hamiltonian. Figure [4] shows 
the spectrum observed for the initial preparation, as well 
as the average of the spectra from the six experiments, 
representing the six terms in /3 of inequality Q with the 
appropriate signs. Fitting the observable spectra, tak- 
ing into consideration the effects of strong coupling, we 
estimate the value of P to be 5.2 ± 0.1, in violation of 
inequality ([2|. The uncertainty on P is propagated from 
the goodness-of-fit figure of merit ascribed to the spectral 
fitting process. 
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FIG. 4: (color online) Summary of experimental results. 
Shown are (in dashed blue) the spectrum produced by the 
initial preparation procedure, Pe ^ establishing a 
reference for e; and (in solid red) the sum of the six spectra 
corresponding to the six terms in (3 of inequality ^ with the 
appropriate signs, scaled by | to compare with the reference. 
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FIG. 5: Numerical simulation results of a simple model of 
decoherence (inset) showing the expected variation of (3 as 
a function of the ratio of the pulse-length, t, to an effective 
dephasing time, T2. The dashed lines indicate bounds on the 
expected performance of the current experiment; for pulse 
length of 1.5ms, and effective decoherence times of 2ms 
T2) and 30ms (^ intrinsic coherence times [H]), the value of 
P is expected to be 1.1 and 5.3 respectively. 

Decoherence, as it is wont to do, causes deviations from 
the idealized closed-system dynamics. To examine its ef- 
fect, we numerically simulate the dynamics of a simple 
model (shown in Figure [5| in which each ideal transfor- 
mation is followed by a symmetric error of a three-fold 
tensor product of a single-qubit dephasing map, A(p), 
given by the operator sum representation p A{p) = 

where Aq = (^l , Ai = J , and 

the parameter r] = 1 — exp {—t/T2) depends on the ratio 
of the pulse- length, to an effective dephasing time, T2. 
Using appropriate estimates [19 of this dephasing time, 
one is able to largely explain the deviation of the experi- 
mental result from the prediction of quantum mechanics 
in ideal conditions. 

Conclusion — We have presented a protocol to di- 
rectly measure correlations between measurement out- 
comes, utilizing an ancillary (probe) two dimensional 
system, with the purpose of testing quantum contextu- 
ality. Conveniently, it can be used directly on ensem- 



bles of quantum systems, without the need for repeated 
projective measurement on single systems. Additionally, 
it can be straightforwardly extended to test similar in- 
equalities on higher-dimensional systems. Our exper- 
imental results demonstrate -under the assumption of 
fair sampling- that a three-qubit deterministic quantum 
computer with one clean qubit reveals correlations be- 
tween measurement outcomes that cannot be explained 
by any NCHV theory. 
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